Light-front quark model analysis of heavy meson radiative decays 
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We present the magnetic dipole(Ml) transitions V — > P7 of various heavy-flavored mesons such 
as (D, D* , D 3 , D*, rj c , J/i/j) and (B , B* , B s , B* using the light-front quark model constrained 

by the variational principle for the QCD-motivated effective Hamiltonian. The weak decay constants 
of heavy mesons and the decay widths for V — * P7 are calculated. The radiative decay for T — > 77^7 
process is found to be very helpful to determine the unmeasured mass of r/b- Our numerical results 
are overall in good agreement with the available experimental data as well as other theoretical model 
calculations. 
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I. INTRODUCTION 

The physics of exclusive heavy meson decays has pro- 
vided very useful testing ground for the precise deter- 
mination of the fundamental parameters of the standard 
model(SM) and the development of a better understand- 
ing of the QCD dynamics. While the experimental tests 
of exclusive heavy meson decays are much easier than 
those of inclusive one, the theoretical understanding of 
exclusive decays is complicated mainly due to the nonper- 
turbative hadronic matrix elements entered in the long 
distance nonperturbative contributions. Since a rigrous 
field-theoretic formulation with a first principle applica- 
tion of QCD to make a reliable estimates of the nonper- 
turbative hadronic matrix elements has not so far been 
possible, most of theoretical efforts have been devoted to 
looking for phenomenological approaches to nonpertur- 
bative QCD dynamics. 

In our previous light-front quark model(LFQM) anal- 
ysis [l[ based on the QCD-motivated effective Hamilto- 
nian, we have analyzed various exlcusive processes such 
as the semileptonic decays between heavy pseudoscalar 
mesons j2| and the rare B — ► K decays [3| and found 
a good agreement with the experimental data. Along 
with those exclusive processes, the magnetic dipole(Ml) 
transitions V(l 3 Si) — > P(l 1 S'o)7 horn the spin-triplet S- 
wave vector (V) mesons to the spin-singlet S'-wave pseu- 
doscalar(P) mesons have also been considered as a valu- 
able testing ground to further constrain the phenomeno- 
logical model of hadrons 0, i, i, 0, 1, 0, M, IS 03, El • 

In this talk we thus investigate the magnetic dipole 
transition among the heavy-flavored mesons such as 
(D,D*,D s ,D*ric, J/ip) and (B,B*,B s ,B* s ,r) b ,T) using 
our LFQM [H, Q ■ Since the experimental data available 
in this heavy-flavored sector are scanty, predictions of a 
model, if found reliable, can be utilized quite fruitfully. In 
addition, we calculate the weak decay constants of heavy 
pseudoscalar and vector mesons. A reliable estimate of 
decay constants is important, as they appear in may pro- 
cesses from which we can extract fundamental quantities 
in the SM such as Cabibbo-Kobayashi-Maskawa matrix 
elements. In our LFQM p], 0, Q, we have implemented 
the variational principle to QCD-motivated effective LF 
Hamiltonian to enable us to analyze the meson mass 
spectra and to find optimized model parameters, which 



are to be used subsequently in the present investigation. 
Such an approach can better constrain the phenomelog- 
ical parameters and establish the extent of applicability 
of our LFQM to wider ranging hadronic phenomena. 

The paper is organized as follows: In Sec. II, we briefly 
describe the formulation of our LFQM 0,0 and the pro- 
cedure of fixing the model parameters using the varia- 
tional principle for the QCD-motivated effective Hamil- 
tonian. The decay constants and radiative V — > P7 de- 
cay widths for heavy-flavored mesons are then uniquely 
determined in our model calculation. In Sec. Ill, the for- 
mulae for the decay constants of pseudoscalar and vector 
mesons as well as the decay widths for V — ► P7 in our 
LFQM are given. To obtain the g 2 -dependent transition 
form factors Fvp(q 2 ) for V — ► P7* transitions, we use 
the Drell-Yan-West q + = frame(i.e. q 2 — — < 0) 
and then analytically continue the spacclike results to 
the timelike q 2 > region by changing to iq± in 
the form factor. The coupling constants gvPj needed 
for the calculations of the decay widths for V — > P7 
can then be determined in the limit as q 2 — > 0, i.e. 
gvp-y — Fyp{q 2 — 0). In Sec. IV, we present our numer- 
ical results and compare with the available experimental 
data as well as other theoretical model predictions. Sum- 
mary and conclusions follow in Sec.V. 



II. MODEL DESCRIPTION 

The key idea in our LFQM [l], [2j for mesons is to treat 
the radial wave function as trial function for the varia- 
tional principle to the QCD-motivated effective Hamil- 
tonian saturating the Fock state expansion by the con- 
stituent quark and antiquark. The QCD-motivated 
Hamiltonian for a description of the ground state meson 
mass spectra is given by 



H q ,\Kt) 



k 2 



= [H 



JJz 



m\ + fc 2 + Vqq 
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(1) 



where k = (k^,fc z ) is the three- momentum of the con- 
stituent quark, M q q is the mass of the meson, and l^^) 
is the meson wave function. In this work, we use two in- 
teraction potentials V q q for the pseudoscalar(0 h ) and 
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vector(l~~) mesons: (1) Coulomb plus harmonic osclla- 
tor(HO), and (2) Coulomb plus linear confining poten- 
tials. In addition, the hyperfine interaction, which is es- 
sential to distinguish vector from pseudoscalar mesons, 
is included for both cases, viz., 



4a. 2 S • S„ 



V„q = Vo + Vhyp = a+Vconf--^ h 77 



3r 



3 m q rriq 



Wcoul, (2) 



where V con f = br(r 2 ) for the linear(HO) potential and 
(Sq-Sg) = l/4(— 3/4) for the vector (pseudoscalar) meson. 

The momentum space light-front wave function of the 
ground state pseudoscalar and vector mesons is given by 



*S(xi,k iX ,Ai) =^ i J J 2 (a; i ,k i x)^(x i ,k i x), (3) 

where ^(x{,kix) is the radial wave function and T^XiX 
is the spin-orbit wave function, which is obtained by the 
interaction independent Melosh transformation from the 
ordinary equal-time static spin-orbit wave function as- 
signed by the quantum numbers J PC . The model wave 
function in Eq. is represented by the Lorentz- invariant 
variables, ir, = pf /P + , kjx = Pix — %iP± and Aj, where 

= (P+,P",Px) = (P° + P 3 ,(M 2 + P2_)/P+,Px) 

is the momentum of the meson M, p% and A, are the 
momenta and the helicities of constituent quarks, respec- 
tively. 

The covariant forms of the spin-orbit wave functions 
for pseudoscalar and vector mesons are given by 
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Ai A2 



e-(pi-P2) 
Ma+m 1 +m 2 



V{P2, A 2 ) 



(4) 



where e M (J z ) is the polarization vector of the vector 



meson [lj], M = M$ - (mi - m 2 ) 2 and Mq is the 
invariant meson mass square M§ defined as Mq = 

^2 k?j_+m " 

= 1 £C, 



The spin-orbit wave functions satisfy 



the following relations E AlA2 ^aYa1^aYa 2 = 1 ' for both 
pseudoscalar and vector mesons. For the radial wave 
function 0, we use the same Gaussian wave function for 
both pseudoscalar and vector mesons 
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^exp(-fc 2 /2/3 2 ), 



(5) 



where (3 is the variational parameter. When the longi- 
tudinal component k z is defined by k z = (x — 1/2) Mq + 
(to 2 — w 2 )/2Mo, the Jacobian of the variable transfor- 
mation {x, kx} — > fc = (kx, fcz) is given by 
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FIG. 1: (Color online). Fit of the ground state meson 
masses[MeV] with the parameters given in Tabled The (p, 7t), 
(77,?/), and (u>, 4>) masses are our input data(red online). 



The normalization factor in Eq. ([5]) is obtained from the 
following normalization of the total wave function, 



dx 



d 2 k 



16tt 3 



^(x,k i± )| 2 = l. 



(7) 



Our variational principle to the QCD-motivated effec- 
tive Hamiltonian first evaluate the expectation value of 
the central Hamiltonian Ho + Vq, i.e. (4>\(Hq + Vb)|0) 
with a trial function <p(xi, kjx) that depends on the vari- 
ational parameters (3 and varies (3 until (4>\(Ho + Vo)\(f>) 
is a minimum. Once these model parameters are fixed, 
then, the mass eigenvalue of each meson is obtained by 
M q q = (4>\(H Q + V q q)\4>). More detailed procedure of de- 
termining the model parameters of light and heavy quark 
sectors can be found in our previous works [l|, [2j . Our 
model parameters (to, (3) for the heavy quark sector ob- 
tained from the linear and HO potential models are sum- 
marized in Table [U 

Our predictions of the ground state meson mass spec- 
tra obtained from the linear and HO potential parameters 
are summarized in Fig.[TJ As one can see, our predictions 
obtained from both linear and H (^parameters are overall 
in good agreement with the data [15| within 6% error. As 
we shall see in our numerical calculations, the radiative 
decay of T — > 77^7 might be useful to determine the mass 
of rib experimentally since the decay width T(T — > r/i,"/) 
is very sensitive to the value of Am(= My — M rjb ) 1 viz. 
Tcx(Ato) 3 . 

The decay constants of pseudoscalar and vector mesons 
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TABLE I: The constituent quark mass[GeV] and the Gaussian paramters /3[GeV] for the linear and HO potentials obtained by 
the variational principle, q = u and d. 



Model 


m q 


m 3 


nic 


mj 




Psc 


Pec 


Pqb Psb 


Pbb 


Linear 


0.22 


0.45 


1.8 


5.2 


0.468 


0.502 


0.651 


0.527 0.571 


1.145 


HO 


0.25 


0.48 


1.8 


5.2 


0.422 


0.469 


0.700 


0.496 0.574 


1.803 



are defined by 



(0\qn 5 q\P) = ifpP*, 
(0\q 7 »q\V(P,h)) = fvM v eP(h), 



(8) 



where the experimental value of vector meson decay con- 
stant fv is extracted from the longitudinal (h = 0) po- 
larization. Using the plus component (/i = +) of the cur- 
rent, one can easily calculate the decay constants and 
the explicit forms of pseudoscalar and vector meson de- 
cay constants are given in [13 ]. 



III. RADIATIVE DECAY WIDTH FOR V -> Pj 

In our LFQM calculation of V — > P7 decay process, we 
shall first analyze the virtual photon(7*) decay process 
so that we calculate the momentum dependent transi- 
tion form factor, F\/p(q 2 ). The lowest-order Feynman 
diagram for V — > P7* process is shown in Fig. [2] where 
the decay from vector meson to pseudoscalar meson and 
virtual photon state is mediated by a quark loop with 
flavors of consituent mass mi and m 2 . 

The transition form factor Fyp(q 2 ) for the magnetic 
dipole decay of vector meson V(P) — ► P(P')y*(q) is de- 
fined as 



The decay form factor Fi/p(q 2 ) can be obtained in the 
q + = frame with the "good" component of currents, 
i.e. /x = +, without encountering zero-mode contribu- 
tions Thus, we shall perform our LFQM calculation 
in the q + = frame, where q 2 = q + q~ — = < 0, 
and then analytically continue the form factor Fyp(q^) 
in the spacelike region to the timelike q 2 > region by 
changing q^ to zqj_ in the form factor. In the calculations 
of the decay form factor Fvp(q 2 ), we use '+ '-component 
of currents and the transverse(/i = ±1) polarization. 

The hadronic matrix element of the plus current, 
(J+) = (P(P')\J+ m \V(P,h = +)) in Eq. © is then ob- 
tained by the convolution formula of the initial and final 
state light-front wave functions: 



(J + ) 



f dx f 

^2 ee j J J£~3 J rf 2 k_L0(s,k'_ L )(?!)(a;,k_ L ; 
xE<I^#T + ^v (10) 



AA 



where k'x = k^ x 2 q_± and eej is the electrical charge for 
j-th quark flavor. Comparing with the right-hand-side of 
Eq. ([!]), i.e. eP + F VP (Q 2 )q R /V2 where q R = q x + iq y , we 
could extract the one-loop integral, I(mi,m 2l q 2 ), which 
is given by 



(P(P')\J^\V(P, h)) = iee^»°e v (P, h)q p P a F VP (q 2 ), 

(9) 

where q = P — P' is the four momentum of the virtual 
photon, e l ,(P,h) is the polarization vector of the initial 
meson with four momentum P and helicity h. The kine- 
matically allowed momentum transfer squared q 2 ranges 
from to 9l 2 nax = (M v 



M P f. 



I{mi,m 2 ,q 2 ) = J J d 2 k 



2 

Mo' 



X1M0K 
(k ± -q ± ) 2 - 



(ii) 



where the primed factors are the functions of final state 
momenta, e.g. M — Mq (a;,k ± ). 

Then, the decay form factor Fyp(q 2 ) is obtained as 




FIG. 2: Lowes-order graph for V — » Pj* transitions. 




F V p(q 2 ) = ei/(mi,TO 2 , q 2 ) + e 2 I(m 2 ,mi,q 2 ). (12) 

The coupling constant gvPj f° r real photon(7) case can 
then be determined in the limit as q 2 — > 0, i.e. gvp-y = 
Fyp(q 2 — 0). The decay width for V — > P7 is given by 



Til' -PD =|<?yp 7 ^, 



(13) 



where a is the fine-structure constant and fc-, 



(M 2 



Mp)/2My is the kinematically allowed energy of the out- 
going photon. 
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FIG. 3: Transition form factors Fvp{q 2 ) for charmed mesons 
radiative decays obtained from the linear parameters. 



IV. NUMERICAL RESULTS 

In our numerical calculations, we use two sets of model 
parameters (m, 0) for the linear and HO confining poten- 
tials given in Table [J to perform, in a way, a parameter- 
free-calculation of decay constants and decay rates for 
heavy pseudoscalar and vector mesons. Although our 
predictions of ground state heavy meson masses are over- 
all in good agreement with the experimental values, we 
use the experimental meson masses except rjb meson in 
the computations of the radiative decay widths to reduce 
possible theoretical uncertainties. Since the rjb mass is 
not measured yet, we use the range Am = M-f — M Vb = 
60 ~ 160 MeV for T->r? 6 7 process [13 ■ 

In Tables [TT] and IIII1 we present our predictions for 
the charmed and bottomed meson decay constants, re- 
spectively, and compare them with other theoretical 
model predictions [ij, QJ, ES E2, El El El EE El 
as well as the experimental data [la, E3, E^, EH 1301 ] . 
Our predictions for the ratios joj Id = 1.18[1.20] and 
fric/fj/ip = 0.91 [0.90] obtained from the linear[HO] pa- 
rameters are in good agreement with the available ex- 
perimental data, (/n s //D)exp. = 1-23 ± 0.11 ± 0.04 [H 
and (/^//j/^cxp. = 0.81 ±0.19 0, E|, respectively. 
Our results for the ratios JbJIb = 1.24[1.32] and 
/b»//b* = 1.23[1.32] obtained from the linear[HO] pa- 
rameters are quite comparable with the recent lattice re- 
sults, 1.20(3)(1) [13 and 1.22(±|) E| for f B *Jfi3* and 
1.17(4)1J H for f B ./f B .. 

We show in Fig. [3] the momentum dependent form fac- 
tors Fvp{q 2 ) for charmed vector meson radiative V — > 
P7* decays obtained from the linear parameters. Since 



the results from the HO parameters are not much differ- 
ent from those of linear ones, we omit them for simplicity. 
The arrows in the figure represent the zero recoil points 
of the final state pseudoscalar meson, i.e. q 2 — g^ ax . We 
have performed the analytical continuation of the decay 
form factors Fyp(q 2 ) from the spacelike region(g 2 < 0) 
to the physical timelike region < q 2 < q 2 nax - The cou- 
pling constant gvp-y at = corresponds to a final 
state pseudoscalar meson recoiling with maximum three- 
momentum in the rest frame of vector meson. The op- 
posite sign of coupling constants for D* + (solid line) and 
D* + (dashed line) decays compared to the charmonium 
J ftp (dot-dashed line) decay indicates that the charmed 
quark contribution is largely destructive in the radiative 
decays of D* + and D* + mesons. The recoil effect, i.e. the 
difference between the zero and the maximum points, is 
not negligible for the D* + — > D + j* decay, while other 
processes may be negligible. The recoil effects for the 
bottomed and bottomonium meson decays are negligible 
due to the very small photon energies. 

In Table [iVl we present our results for the decay 
widths and branching ratios together with the available 
experimental data. The errors in our results for the 
decay widths and branching ratios come from the un- 
certainties of the experimental mass values and exper- 
imental mass values plus the full widths, respectively. 
Our results of the branching ratios Br(J/%p — ► 7^7) = 
1.80 ± 0.10[1.76 ± 0.10]% and Br(L>* -> D+j) = 0.93 ± 
0.31[1.00±0.34]% obtained from the linear[HO] parame- 
ters are in agreement with the experimental data [Til ]. 
Br(J/i/> -► r?c7)ox P = (1.3 ± 0.4)% and Br(D* -> 
-D + 7) cxp = (1.6 ± 0.4)% within the error bars. For 
the T — * T7&7 process, our predictions for the decay 
width and branching ratio obtained from the linear [HO] 
parameters are T(T -> rj b j) = 45l3g[42l||] eV and 

Br(T -> rj b j) = (8.4l£ , j 8 )[7.7l& ] x 10~ 4 , where the 
lower, central, and upper values correspond to Am = 60 
MeV, 110 MeV, and 160 MeV, respectively. The decay 
width T(T — > r/b-f) is found to be very sensitive to Am 
because it is proportional to (Am) 3 . Other model cal- 
culations for the T(15) radiative Ml decay rates can be 
found in Ref. [H|. 

In Fig. we show the dependence of T(T — + 77^7) 
on Am compared with other theoretical model calcula- 
tions [13] • As one can see from Fig. QJ our prediction for 
the dependence of T(T — * 77^7) on Am is quite consitent 
with other theoretical predictions for various Am (32j . 



V. SUMMARY AND DISCUSSION 

In this work, we investigated the weak decay con- 
stants and the magnetic dipole V — ► Pj decays of 
heavy-flavored mesons such as (D, D* , D s , D*,rj c , J/tp) 
and (B, B* , B s , B*,rjb,T) using the LFQM constrained 
by the variational principle for the QCD-motivated effec- 
tive Hamiltonian. Our model parameters obtained from 
the variational principle uniquely determine the above 
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TABLE II: Charmed meson decay constants(in unit of MeV) obtained from the linear[HO] parameters. 





So 


Id* 


Sn s 


hi 


he 


Sj/i> 


Linear [HO] 


211(194] 


254 [228] 


248(233] 


290 [268] 


326 [354] 


360(395] 


ij&ttice |io] 


Oil J_ 1 A^~' 2 
Zll It l i ±_22 


Z40 zt ^U_2 


OQ1 _|_ 1 O+S 
ZOl zt 








QCD [19] 


201±3±17 




249 ± 3 ± 16 








Sum-rules [20] 


204 ± 20 




235 ± 24 








BS [21] 


230 ± 25 


340 ± 23 


248 ± 27 


375 ± 24 


292 ± 25 


459 ± 28 


QM [22] 


240 ± 20 




290 ± 20 








RQM [23] 


234 


310 


268 


315 






Exp. 


222.6 ± I6.7+3 4 f 27 l 




274 ± 13 ± 7 [28] 




335 ± 75 [29] 


416 ±6 [15] 



TABLE III: Bottomed meson decay constants(in unit of MeV) obtained from the linear[HO] parameters. 





Sb 


Sb* 


Sb. 


Sb: 


Svb 


fr 


Linear [HO] 


189(180] 


204 [193] 


234 [237] 


250 [254] 


507[897] 


529(983] 


Lattice [18] 


179 ± 18+f 


196 ± 24+if 


204 ± 16+.jf 


229 ± 20+ 4 6 






QCD [24] 


216 ±22 




259 ± 32 








[25] 


189 ± 27 




230 ± 30 








Sum-rules [26] 


210 ± 19 




244 ± 21 








[20] 


203 ± 23 




236 ± 30 








BS [21] 


196 ± 29 


238 ± 18 


216 ±32 


272 ± 20 




498 ± 20 


QM [22] 


155 ± 15 




210 ± 20 








RQM [23] 


189 


219 


218 


251 






Exp. 


229+^+37 [3Q] 










715 ±5 [15] 
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FIG. 4: The dependence of F(T -> 7767) on Am = Mr - M, 
compared with other theoretical model calculations [32]. 



good agreement with the available experimental data as 
well as other theoretical model calculations. Our nu- 
merical results of the decay widths for J/ip — > r\ c ^ and 
D* + — ► Z3+7 fall within the experimental error bars. 
Our predictions for the branching ratios for the bot- 
tomed and bottomed-strange mesons are quite compar- 
ible with other theoretical model predictions. For the 
radiative decay of the bottomonium, we find that the de- 
cay widths r(T — > rjb'y) is very sensitive to the value of 
ArTi = My — M r/b . This sensitivity for the bottomonium 
radiative decay may help to determine the mass of r\b 
exp er iment ally. 

Since the form fator Fvp(q 2 ) of vector meson radia- 
tive decay V — * Pj* presented in this work is precisely 
analogous to the vector current form factor g(q 2 ) in weak 
decay of ground state pseudoscalar meson to ground state 
vector meson, the ability of our model to describe such 
decay is therefore relevant to the reliability of the model 
for the weak decay. Consideration on such exclusive weak 
decays in our LFQM is underway. Although our previous 
LFQM [H,[| and this analyses did not include the heavy 
mesons comprising both c and b quarks such as B c and 
B* , the extension of our LFQM to these mesons will be 
explored in our future communication. 



nonperturbative quantities. This approach can establish 
the extent of applicability of our LFQM to wider ranging 
hadronic phenomena. 

Our predictions of mass spectra and decay constants 
for heavy pseudoscalar and vector mesons are overall in 
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TABLE IV: Decay widths and branching ratios for radiative V — > P7 decays obtained from our linear[HO] model parameters. 
We used M Vb = 9353 ± 50 MeV for T -> 17,, 7 decay. 



Decay mode 


1 [Kevj 




r>r cxp |lo| 


J/ip — * rj c "i 


1.69 ± 0.05(1.65 ±0.05] 


(1.80 ±0.10)[1.76± 0.10]% 


(1.3 ±0.4)% 


D *+ _^ D+1 


0.90 ± 0.02(0.96 ±0.02] 


(0.93 ±0.31)[1.00± 0.34]% 


(1.6 ±0.4)% 


D* u -> _D U 7 


20.0 ±0.3(21.0 ±0.3] 




(38.1 ±2.9)% 


or -» o s + 7 


0.18 ±0.01(0.17 ±0.01] 




(94.2 ± 0.7)% 




0.40 ± 0.03(0.40 ± 0.03] 






B* u -» B u 7 


0.13 ±0.01(0.13 ±0.01] 






5.*° - B s u 7 


0.068 ±0.017(0.064 ±0.016] 






T — > ^7 


0.045i u ;^[0.042tr Q f 6 ] 


(8.4i£ b )[7.7i^] x 10- 4 
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